arXiv: 1501.07365v2 [cs.RO] 29Jun2015 


7R Darboux Linkages by Factorization of Motion Polynomials 


Z. Li 
RICAM 

Austrian Academy of Sciences 
Linz, Austria 


J. Schicho 
RISC 

Johannes Kepler University 
Linz, Austria 

June 30, 2015 


H.-R Schrocker 
Unit Geometry and CAD 
University of Innsbruck 
Innsbruck, Austria 


Abstract 

In this paper, we construct two types of 7R closed single loop linkages by combining different factorizations of a general 
(non-vertical) Darboux motion. These factorizations are obtained by extensions of a factorization algorithm for a generic 
rational motion. The first type of 7R linkages has several one-dimensional configuration components and one of them 
corresponds to the Darboux motion. The other type is a 7R linkage with two degrees of freedom and without one-dimensional 
component. The Darboux motion is a curve in an irreducible two dimensional configuration component. 
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1 Introduction 

In 1881, G. Darboux determined all one-parametric spatial 
motions with only planar curves as point paths [1]. These 
properties determine a unique motion which is nowadays 
named after Darboux. It is the composition of a planar 
elliptic motion with a suitably oscillating translation perpen¬ 
dicular to that plane. If the elliptic motion degenerates to 
a rotation, the Darboux motion is called vertical, if this is 
not the case, we call the Darboux motion non-vertical or 
general. In both cases, the generic point paths are ellipses. 
The Darboux motion is a special SchoenflieB motion and we 
call the direction perpendicular to the planar elliptic motion 
component the axis direction. These and other properties are 
well known in literature, see for example [2] or [3, Chapter 9, 
§3]. 

In this article, we construct 7R linkages with one degree 
of freedom such that one link during (one component of) the 
motion performs a general Darboux motion. These results 
complement recently presented linkages that generate ver¬ 
tical Darboux motions [4]. The beginnings of this research 
can be found in [5] where we investigated spatial linkages 
with a straight line trajectory. The linkages we present in this 
article are novel. Unlike the linkages of [4], we generate a 
general (non-vertical) Darboux motion and our linkages are 
free of prismatic or cylindrical joints. One of our linkages 
has already be found in [5]. Our description here is more 
general and detailed. 

Our linkage construction is based on the factorization 
of rational motions which has been first introduced in [6]. 
In the dual quaternion model of rigid body displacements, 


the Darboux motion is parameterized by a cubic motion 
polynomial. In general, a cubic motion can be generated 
by 6R linkages [7-9] but specialties of the Darboux mo¬ 
tion prevent application of the factorization algorithm for 
generic motion polynomials. Therefore, we resort to new 
factorization techniques for non-generic motion polynomi¬ 
als [10-12], i.e., we multiply the motion polynomial with a 
real polynomial which divides the primal part of the motion 
polynomial. The new motion polynomial can be factored 
in different ways which allows us to obtain new closed 7R 
linkages by combining different factorizations. A unified al¬ 
gorithm for factorization of generic and non-generic motion 
polynomials is introduced in [11]. Here, we just mention the 
basic idea of this algorithm: By factoring out linear motion 
polynomial factors (from the left and from the right) we 
create a generic motion polynomial that can be factored by 
the original algorithm of [6]. 

We continue this text with a quick introduction to the dual 
quaternion model of SE(3), the group of rigid body displace¬ 
ments, and to motion polynomials and their factorization 
(Section 2). In Section 3 we describe the Darboux motion in 
the dual quaternion model, present different factorizations, 
and combine them to form 7R linkages. 

2 Preliminaries 

The investigations in this article are based on the dual quater¬ 
nion model of the group SE(3) of rigid body displacements. 
We assume familiarity with dual quaternions and their rela¬ 
tion to kinematics but we provide a quick introduction of 
concepts that cannot be found in standard textbooks. For 
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more information on basics of dual quaternions and kine¬ 
matics we refer to [13-15]. A short presentation of the topic 
with close relations to this article can also be found in [16]. 

2.1 Dual quaternions 

A dual quaternion is an object of the shape 

h = ho + h\i + h 2 \ + / 13 k + £(7*4 + h$\ + + / 17 k) ( 1 ) 

with real numbers The associative but non- 

commutative multiplication of dual quaternions is defined 
via the relations 

i 2 =j 2 =k 2 =ijk=-l, £ 2 =0, 

\e = ei, j£ = £j, k£ = £k. 

For example 

(1 -ei)(j + ek) =j + £k-eij-£ 2 ik 
= j + £k-ek-0 = j. 

Conjugate dual quaternion and dual quaternion norm are 

h = ho-hii-h 2 i-h 3 k + £(h 4 — h$i — — h 2 k), 

\\h\\=hh. 

The primal part of h is p = ho + h\\ + h 2 ] + h 3 k and the dual 
part is d = h 4 + h$\ + + h 2 k. Primal and dual part are 

(ordinary) quaternions. The set of quaternions is denoted by 
H, the set of dual quaternions by BH := B (Sir H. Writing 
h = p + £d with p,d E H, we have \\h\\ = pp J r£(pq J rqp). 
Note that both pp and pq + qp are real numbers and the 
norm itself is a dual number. 

If the dual quaternion h = p-\-£d has a real norm (pq + 
qp = 0), its action on a point [xo,x\,x 2 ,x 3 \ of real projective 
three-space is defined via 

= (2) 
PP 

where x = xo + x\i + x 2 } + x 3 k, y = y 0 + y ii + yii + ^k and 
the image point has projective coordinates [yoAiA2A3]- 
The map (2) describes a rigid body displacement and the 
composition of such displacements corresponds to dual 
quaternion multiplication. Note that h and Xh with A G 
M \ {0} describe the same map. 

The dual quaternion h as in (1) describes a rotation or 
translation, if and only if h 4 = 0. Among these, transla¬ 
tions are characterized by h\ = h 2 = h 3 = 0. We speak 
accordingly of rotation or translation quaternions. In 
the former case, the revolute axis has Plucker coordinates 
[/ii,/* 2 ,/* 3 , —^ 5 , —^ 6 , —hj\. In particular, the axis direction 
is (h\,h 2 M)- 


2.2 Motion polynomials 

Now we study polynomials 

C = c n t n + • • • -\-c\t + co (5) 

in the indeterminate t and with dual quaternion coefficients 
c n ,... , co G ML The set of these quaternions is denoted by 
Multiplication of polynomials in MI [t] is defined 
by the convention that indeterminate t and coefficients com¬ 
mute. The conjugate polynomial C is the polynomial with 
conjugate coefficients c^,..., cfi. 

A polynomial C G MI|Y] is called a motion polynomial, 
if its leading coefficient c n is invertible (has non-zero pri¬ 
mal part) and the so-called norm polynomial CC has real 
coefficients. In this case, it parameterizes a rational motion 
via the map (2). Note that C and QC parameterize the same 
rational motion if Q is a real polynomial without zeros. This 
observation will be important in the next section. 

The linear polynomial C = t — h with h as in (1) is a mo¬ 
tion polynomial if and only if J 14 = h\h 3 + h 2 h^ + h 3 h 2 
0 . Assuming (h\,h 2 ,h 3 ) 7 ^ (0,0,0), the displacements de¬ 
scribed by C for varying t are rotations about a fixed axis but 
with different revolute angle. The axis’ Plucker coordinates 
are [h\,h 2 ,h 3 ,—h 3 ,—h 6 ,—h 2 \ and (h\,h 2 ,h 3 ) is a direction 
vector of the axis. 

A generic motion polynomial C of degree n admits n\ 
factorizations of the shape 

C=(t-hi)~-(t-h n ) (4) 

with rotation quaternions h\,...,h n . The precise statement 
and a factorization algorithm can be found in [ 6 ]. Here is 
suffices to say that a Darboux motion is not generic and the 
general theory of [ 6 ] gives no information about existence 
or non-existence of factorizations. 

Equation (4) admits an important kinematic interpretation: 
Each factorization corresponds to a one-parametric motion 
of an open chain of n revolute joints whose end-effector fol¬ 
lows the motion parameterized by C. The revolute axes are 
determined by the rotation quaternions h\ y . , . 1 h n . Different 
factorizations may be combined to form closed loop linkages 
with this coupler motion. 

In the remainder of this text we need auxiliary results 
related to the factorization of motion polynomials. Their 
proofs can be found in [ 6 ]: 

Lemma 1 (Polynomial division). If C and D are polyno¬ 
mials in BH^] and the leading coefficient of D is 1, there 
exist polynomials Q,R G BH [t] such that C = QD -\-R and 
deg R < degD. 

Note that the computation of quotient Q and remainder R is 
possible by straightforward polynomial long division in a 
non-commutative setting. 
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Lemma 2 (Zeros and right factors). The dual quaternion 
h is a zero of the polynomial C E BH [t] if and only if there 
exists a polynomial Q E BIHI|Y] such that C = Q(t — h). 

The statement of Lemma 2 is less trivial than the corre¬ 
sponding statement for real polynomials. In particular, with 
C as in (3), the value of C(h) (and hence also the term “zero 

of C”) is defined as c n h n H-j- c\h + co- It is not admissible 

to simply plug h in a factorized representation of the shape 
(4). With C = (t — k)(t — h) and h,k E DM [t] we have, for 
example, C = t 2 — (h + k)t+ kh and thus 

C(h) = -kh + kh = 0 but C(k) =hk + kh^0. 

Lemma 2 asserts that the “usual” relation between zeros 
linear factors holds at least for right factors. 

3 Factorizations of the Darboux motion 


In [5], the authors gave two types of factorizations for the 
motion C in Equation (6). The first type is with three linear 
motion polynomials (FI), the second type (FII) is with five 
linear motion polynomials whose product does not equal 
C but C times a real quadratic polynomial P. This makes 
an algebraic difference but does not change the kinematics: 
Both C and PC parameterize the same motion. Moreover, 
two successive linear factors in FII are identical. Hence, 
it gives rise to an open 4R chain which can generate the 
Darboux motion. In combination with the 3R chain obtained 
from FI we obtain a 7R linkage. 

In this section, we recall the construction of [5] and we 
give two further factorizations (Fill and FIV) of (6). Similar 
to FII, these factorizations also use five linear motion polyno¬ 
mials but give only four axes because two successive linear 
motion polynomials are same. Combining it with the already 
known 3R chain, we obtain further 7R Darboux linkages. 


Following [3, Equation (3.4)], the general Darboux motion 
is given by the parametric equations 

X = vcos<p — ysintp 
Y = vsintp +ycos<jP -\-asin(p 
Z = z + h sin (p + c( 1 — cos <p). 


3.1 First factorization FI 

We want to find linear motion polynomials Q\, Q 2 , Q 3 such 
that C = <2i 22*23 • By Lemma 2, a necessary and sufficient 
condition for 23 = t — # 3,^3 E BH, to be a right factor of C 
is C(qf) = 0. A straightforward calculation shows that qs is 
necessarily of the shape 


Here, (p E [— 71 , 7 t) is the motion parameter, a , h , c are real 
constants and x, y, z resp. X , Y, Z are coordinates in the 
moving and the fixed frame. For a = 0, the motion is a 
vertical Darboux motion. In this paper we exclude this case 
and assume a f 0. 

Using the conversion formulas between matrix and dual 
quaternion representation of rigid body displacements (see 
for example [15]), we find the polynomial parameterization 

Co = (k + ce)t 3 + (1 — e(ai — ck + b))t 2 

+ (k-e(flj + fck))f + l (5) 


<73 = k + £(vi + wj), v,weM. 

Next, we consider the quotient polynomial Q = P + eD , 
defined by C = 223- Us primal part is P = t 2 + 1. Thus, the 
motion parameterized by 2 is a curvilinear translation. Since 
2 is quadratic and P has no real zero, the trajectories are 
congruent ellipses. This motion admits a factorization Q = 
2122 of the required shape only if it is a circular translation 
(see [16] for a detailed proof of this statement). This puts a 
condition on # 3 , namely 

v=—a~ 1 bc, w = — (la)~ l (a 2 -\-b 2 — c 2 ). (7) 


for the Darboux motion. Left-dividing this motion polyno¬ 
mial by the leading coefficient k + c£ (this amounts to a 
coordinate change in the fixed frame), we arrive at 

C = t 3 — (k — e(aj — bk))t 2 -{- (1 — e(b-\-ai — ck))t — k + c£. 

( 6 ) 

Our further investigations will be based on this parametric 
representation of the Darboux motion. 

Note that the primal part of C is (t 2 + l)(t — k). The 
presence of a real factor violates the assumptions of [6, The¬ 
orem 1] (existence of factorizations) and prevents application 
of Algorithm 1 (computation of factorizations) of that paper. 
Nonetheless, we will see below that C admits infinitely many 
factorizations with three linear motion polynomials. More¬ 
over, suitable polynomial multiples of C admit factorizations 
as well and parameterize the same Darboux motion. This 
we will use for our construction of 7R Darboux linkages. 


If (7) is satisfied, the remainder Q is a circular translation 
and can be factored in infinitely many ways as Q = Q 1 Q 2 
with linear rotation polynomials 2i> 22- Each pair of these 
factorizations correspond to a parallelogram linkage that 
generates the circular translation. One example of a factor¬ 
ization of C is 


be . a 2 -f c 2 — b 2 . 

2i =t + E - i£H --- je-feke, 

a 2 a 

Q 2 —t E, 

be . a 2 + b 2 — c 2 . 

23 —t — k H-1£ H---J£ 

a 2 a 


( 8 ) 


where 


E = 


lac 


lab 


i 2 + b 2 + c 2 a 2 + b 2 + c 2J a 2 j rb 2 j rc 2 


a 2 — b 2 — c 2 

J + .9 , r,9 , _9 k - (9) 


In particular, we have shown 
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Proposition 3. The general Darhoux motion parameterized 
by ( 6 ) can be decomposed into a rotation about an axis par¬ 
allel to the Darboux motion’s axis direction and a circular 
translation in a plane orthogonal to the vector quaternion D 
of (9). 

Combining two of the infinitely many factorizations with 
three linear factors gives a parallelogram linkage to whose 
coupler a dangling link (corresponding to the right factor Q 3 
which is the same for all factorizations) is attached. This 
6 R linkage has two trivial degrees of freedom and is not 
particularly interesting. 

3.2 Second factorization FII 

Now we construct a further factorization FII of the general 
Darboux motion. We do, however, alter the parameteri¬ 
zation ( 6 ) and multiply it with the quadratic polynomial 
P := t 2 + 1. This changes the motion polynomial but not the 
parameterized motion and gives additional freedom to find 
factorizations. In contrast to FI, we have more than three 
linear factors. 

We begin by setting the right factor to Q 4 := t — k and 
then use polynomial division (Lemma 1) to compute 

Ci = t 2 + e(ai — bk)t + 1 + ce k. 

such that C = C 1 Q 4 . The polynomial PC\ can be written 
as the product of four linear motion polynomials. Again, 
there exist infinitely many factorizations. In order to keep 
the number of joints small, we choose one with identical 
middle factors, i.e., PC\ = QiQ\Q$, where 

^ . a+c . b , 

Qi — t + H —— — 2^’ 

Qe=t-i, 

a—c . b . 

Q 5 =t + i+-^-£)--£k. 

The multiplicity of the middle factor Q 5 allows us to make 
a 3R chain for the motion parameterized by PC \. Together 
with C = Q 1 Q 2 Q 3 and Q\, we can get a 7R linkage. It can 
be seen that the axes of Q 1 , Q 2 are parallel, as are the axes 
of 23 , 24 and 25 , 26 , 27 - 

For a concrete example, a Grobner basis computation re¬ 
veals that the configuration set of this 7R linkage is one 
dimensional [17]. Note that its configuration curve con¬ 
sists of different components and only one corresponds to 
the Darboux motion parameterized by C. We also want to 
mention that the construction of the second factorization 
FII leaves many degrees of freedom. We just presented one 
factorization of this type. 

3.3 Third factorization Fill 

In order to find a third factorization Fill that is suitable for 
linkage construction, we proceed in a similar fashion to case 


FII, multiplying the motion polynomial C of Equation (6) 
with P = t 2 + 1. But we do so prior to splitting off a right 
factor The polynomial PC has infinitely many factorizations 
into products of five linear motion polynomials, one of them 
being of the shape CP = Q'-jQ'^Qf ■ 

We start by finding a suitable factor Q' 4 . Because Q' 4 Q ' 4 is 
a factor of (CP) (CP) = P 5 , we have the necessary constraint 
Q\Q , 4 = P but still have infinitely many choices. We take 

24 = t — k — xie — yj£ 

with x,y E M. After division of PC by (Q 4 ) 2 (Lemma 1), 
we are left with the cubic motion polynomial C 2 defined by 
CP = C 2 Q 4 , i.e., 


C 2 = t 3 + (k + e(aj — bk + 2 xi + 2 yf))t 2 

+ (1 + e(ai + ck + 2yi — 2xj + b))t + k — c£. 


This new cubic motion polynomial C 2 parameterizes again 
a Darboux motion but the parameterization is not in the 
standard form (6). Nonetheless, there exist again infinitely 
many factorizations similar to (8). For example, we have 

C2 = QiQ'eQ's where 


Qs — t T - k + ^ 


b 2 x — abc — 2 bey — c 2 x 


T 2 + Ax 2 
/ ab 2 — ac 2 + 2 b 2 y + 4 bex — 2 c 2 y 


+x ie 


V 


2(r 2 +4x 2 ) 


a 

+ - 


+y)j£, 

, 2 (ac — 2 bx-\- 2 cy) . 2 (ab-\- 2 by -\- 2 cx) . 

° 6 = t ~ £ 2 + c 2 + T 2 + 4x 21 b 2 c 2 T 2 4x 2 J 


+ 


T 2 — b 2 — c 2 + 4x 2 
b 2 + c 2 + T 2 + 4x 2 


K 


Qi = q: 


—t Sx + bcT + 4 x(ay + x 2 + y 2 ). 


r 2 + 4x 2 


\£ 


(S-\-4x 2 )aT -\-4(ay(a(3y 4-a) +2y 2 ) —abex), 
3 4 77^5 i ^ o \ 


2a(T 2 + 4x 2 ) 


-feke, 


and we abbreviated S = a 2 — b 2 -\-c 2 and T = a + 2y. 

Because C = 2 i 2223 and PC = C 2 Q 4 = Q' 7 Q' 6 Q' 5 Q 4 
parameterize the same motion, we can combine these two 
factorizations to form a 7R linkage where each rotation is 
defined by 2y, 2^ 2s, 2i, 23, 22, 2i- It can be seen that 
the axes of 2i, Qi are parallel, as are the axes of 23, 2^, 25 
and 2v> 26- 

As revealed by a Grobner basis computation, the configu¬ 
ration space of a concrete numeric example of this 7R link¬ 
age is really a curve [17]. Thus, we have indeed constructed 
another 7R linkage whose coupler motion is a non-vertical 
Darboux motion. Note that the configuration curve contains 
several components, not all of them rational. One component 
corresponds to the rational curve parameterized by C. 

In Figure 1, we present nine configurations (the first one 
and the last one are from same configuration) of this linkage 
in an orthographic projection parallel to k. 
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Figure 1: The first 7R linkage generates a non-vertical Darboux motion. 
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3.4 Fourth factorization FIV 

Finally, we present an exceptional case of the 7R linkage ob¬ 
tained from the third factorization Fill. Among the infinitely 
many choices for a, b, c, x, y we choose the special case 
a = 1, b = 2, c = 0, x = 0, y = 0. We then have 

Q'l = t + fj - - 2k£, 

Q" = t — k + -j£, 

G4=*-k> 
fi5=^ + k+ 2^’ 

2 " = t ~ - ^> £ - 2k£ - 

We can combine this factorization to form a 7R linkage 
where each rotation is defined by Q", Q'l, Q'l, Q'l, Q 3 , 2^ 
2'/. It can be seen that the axes of Q'{, Q'l, Q'l, Q'l are 
parallel, as are the axes of Q'l, Q'l, Q'l. Since four neighbor¬ 
ing axes are parallel, this linkages contains a planar 4-bar 
linkage. Furthermore, there are two pairs of parallel 3R 
sub chains which give us two Sarrus linkages if we fix the 
remaining joint. Hence, this 7R linkage has at least two 
degrees of freedom. 

A Grobner basis computation confirms that the configu¬ 
ration space is indeed two-dimensional [17]. Moreover, a 
decomposition of the configuration variety shows that the 
linkage is not kinematotropic, which means that there are 
components of different dimension [18]: It consists of three 
irreducible 2-dimensional components. Two components are 
rational. We do not know whether the other two-dimensional 
component which contains the rational curve parameterized 
by C is rational or not. 

In Figure 2, we present eleven configurations (the first 
one and the last one are the same) of this linkage in an or¬ 
thographic projection parallel to the direction of k. Because 
the axes of Q'[, Q'l, Q'l, and Q'l are parallel, the remaining 
axes are depicted as points throughout the motion. We can 
observe the parallelity of axes and possible singularities (in¬ 
tersections). The two configurations where the axes of Q'l 
and Q'l coincide give us two rational configuration sets of 
dimension two: A planar 4R loop (the coupler motion is a 
circular translation) and a single rotation rational component. 
None of these two contains the Darboux motion. 

4 Conclusions 

Using factorization of motion polynomials in non-generic 
cases, we constructed several examples of new 7R Darboux 


linkages. The main idea is to find an irreducible real poly¬ 
nomial P and factor PC instead of C. The obtained 7R 
linkages generate the general (non-vertical) Darboux motion 
and exhibit some interesting specialities like parallel axes 
or circular translations as relative motions between certain 
links. One can replace two R-joints by other R-joints which 
generate same circular translations. 

Verifying the validity of the presented factorizations and 
corresponding linkages is trivial but, admittedly, we did not 
explain in detail how to actually find the factorizations. The 
factorization theory for non-generic motion polynomials, 
including theoretical results, algorithms and upper bounds 
for the number of factors, is currently being worked out. Its 
presentation is left to future publications. 
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